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Abstract. This work presents a multiple scattering formulation of two dimensional 
acoustic metamaterials. It is shown that in the low frequency limit multiple scattering 
allows us to define frequency-dependent effective acoustic parameters for arrays of 
both ordered and disordered cylinders. This formulation can lead to both positive 
and negative acoustic parameters, where the acoustic parameters are the scalar bulk 
modulus and the tensorial mass density and, therefore, anisotropic wave propagation 
is allowed with both positive or negative refraction index. It is also shown that the 
surface fields on the scatterer are the main responsible of the anomalous behavior 
of the effective medium, therefore complex scatterers can be used to engineer the 
frequency response of the effective medium, and some examples of application to 
different scatterers are given. Finally, the theory is extended to electromagnetic wave 
propagation, where Mie resonances are found to be the responsible of the metamaterial 
behavior. As an application, it is shown that it is possible to obtain metamaterials 
with negative permeability and permittivity tensors by arrays of all-dielectric cylinders 
and that anisotropic cylinders can tune the frequency response of these resonances. 
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1. Introduction 

Fluid-like metamaterials or metafluids with negative constitutive parameters offer a new 
insight into acoustic wave propagation. Single negative metamaterials (SNM), in which 
either the mass density or the bulk modulus is negative [Ill2l[3], can be used, for example, 
designing of surface-like acoustic lens to overcome the diffraction limit [H |5] or for the 
design of acoustic panels [6]. Double negative metamaterials (DNM)[71 [U [9] present 
negative refraction [101 [HI 121 US] and, as it is well known from electromagnetic wave 
theory, they can also be used to increase the resolution of conventional lens p3| [T5| IT6] . 
In general, anisotropic fluid-like metamaterials with acoustic parameters both positive 
and negative are necessary in the field of transformation acoustics for the design of 
several types of acoustic devices based on sound propagation [TTl [HI [191 (SHI [21] . 

The existence of frequency ranges in which the effective medium presents negative 
constitutive parameters is related with the subwavelength resonances of the individual 
scatterers that constitute the metamaterial, being these due to soft-scatterer resonances 
[71 [22l [23] or due to Helmholtz-like resonances [231 [23 [211 [H [2] . The same phenomenon 
is found in electromagnetic waves under the name of Mie resonances [271 [2H1 [29l [30], and 
they present an alternate way of design electromagnetic metamaterials to that offered 
by split ring resonators [31] or metallodielectric composites [13], which have been the 
dominant structures so far. Therefore, metamaterials based on the resonances of the 
individual scatterers are important not only for acoustic but also for electromagnetic 
metamaterials. 

It is known that a monopolar resonance in the scatterer is the responsible of the 
negative bulk modulus, and that a dipolar one is the responsible of the negative mass 
density [7]. However, the full effect of the ensemble of scatterers that constitute the 
effective medium has been partially explained only, as multiple scattering effects or 
anisotropic lattices have not been considered yet. 

In this work, we give a comprehensive description of multiple of acoustic 
metamaterials by using a multiple scattering formulation. It is based on our previous 
results on homogenization of sonic crystals [321 [SSI [M] . This formulation describes, in 
the long wavelength limit, an ensemble of orderer or disordered scatterers as an effective 
medium with frequency-dependent acoustic parameters, which are shown to be negative 
in certain frequency regions. The frequency-dependent parameters are given in terms of 
both the lattice symmetry and the surface fields on the scatterers, showing that these 
fields are the quantities that we have to manage in order to engineer the frequency 
response of our effective medium. 

Therefore this formulation allow us, from one side, summarize all the previous 
results regarding SNM and DNM, from the other side, extend the theory to any 
type of radially symmetric scatterer and to non symmetric lattices, like rectangular 
arrangements in place of square or hexagonal ones, which are the more usually studied. 
Also, we show that the regions in which the metamaterial presents divergent or negative 
acoustic parameters is mainly a function of the surface fields of the scatterers, opening 
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therefore the possibihty of engineer the scatterer in order to increase the frequency 
region in which the anomalous behavior occurs. 

The paper is organized as follows: In section [2]the concept of quasi-static resonance 
is introduced, showing how these resonances lead to scatterers with locally negative 
parameters. Section Appendix A analyzes several examples of scatterers with locally 
negative parameters, showing that complex scatterers like fluid-like shells or anisotropic 
fluid-like cylinders can tune the metamaterial behavior of the scatterer. After that, in 
Section H] the multiple scattering formulation is presented, analyzing first the case in 
which the scatterers have small radius (low filling fraction) and, later, the more general 
case, showing how anisotropy appears. In this section the theory is extended to include 
multipolar effects, but we see that they are not very important in principle. Finally, the 
application of the theory to electromagnetic waves is explained in section [5l The paper 
ends with a summary section. 



2. Quasi-Static Resonances and Locally Negative Parameters 

Homogenization theories for assembles of scatterers are based on the low wavenumber 
(long wavelength) expansions of the fields in both the background and the scatterers. 
When working with metamaterials we assume that the wavenumber in the background 
is asymptotically small though we let the wavenumber inside the scatterer still be finite. 
Physically it means that outside the scatterers the wave field propagates through an 
effective medium but it is still allowed to the scatterers to have complex internal 
scattering processes, which will lead to locally (i.e., in a narrow frequency region) 
negative parameters, as will be explained in the following sections. 

The simpler example of these scatterers, and the first studied in the next section, 
is the homogeneous fiuid-like scatterer. If the speed of sound inside this scatterer is 
much smaller than that of the background, Cq << Cb, we will have that, for a given 
frequency u, the wavelength inside the scatterer Xa will also be much smaller than that 
in the background << Xb- Thus, outside the scatterer the field will be a function of 
kb = u/cb, which is a slow oscillating function, while inside the scatterer the fields will 
be a function of ka = ^^/ca and, therefore, it is a rapidly oscillating function. As we 
are in the low frequency limit we expect the medium behave as a homogeneous effective 
medium with constant parameters, but, due to the fields inside the scatterer, we will 
find that our effective medium has frequency- dependent parameters. 

Next section analyzes this effect rigorously and for several type of scatterers, 
showing how complex scatterers present a frequency behavior that, as will be seen 
in Section HJ can lead to metamaterials with negative constitutive parameters. 
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3. Acoustic Scatterers with Locally Negative Parameters 

The wave equation for the pressure field in an inhomogeneous fiuid is given by [35] 



V(p-^(r)VP(r)) + -^P(r) = (1) 
Bir) 



where p{r) and B{r) are the mass density and bulk modulus, respectively, of the fiuid 
material. The problem considered here is reduced to points in the x — y plane, that in 
polar coordinates are r = {r,6). We also assume that a radially symmetric scatterer of 
radius Ra with some inhomogeneous parameters p(r) and B{r) is embedded into a fiuid 
background with acoustic parameters pi, and Bf,. 

This is a canonical scattering problem whose solution outside the scatterer is given 
in terms of Bessel and Hankel functions [35] , 

P{r,e-u) = Y,Al[Jg{hr)+TgHg{hr)]e"^' , r > R, (2) 
q 

with kl = u'^pb/Bb. The coefficients are determined by the incident field, and 
the response of the scatterer is described by the matrix elements Tq. This matrix is 
obtained by solving the wave equation ([1]) inside the scatterer and applying boundary 
conditions at r = Ra, which are the continuity of the pressure field and that of the 
normal component of the particle velocity, 

PiRt) = PiR-) (3a) 

-drP{Rt) = -^drP{R-). {3b) 

Pb P[Ra) 

Since the scatterer is radially symmetric, and the parameters p and B depend only 
on the radial coordinate, the field inside the scatterer is expressed in a Fourier series of 
the form 

P{T,e-u) = Y,B,{uj)^,{r-u)e^', (4) 

where the eigenfunctions ■ipq{r]u) are solutions of the radial part of in cylindrical 
coordinates, 

P^d. {j^frMr; -)) + (-^||^ - ^) Ur; = 0. (5) 

From this equation, after applying the boundary conditions, the general form for the T 
matrix is obtained 

X.H'^ihRa) - HqihRa) ' P, dr^Ra, ^) ' ^' 

This T matrix allows to distinguish the contribution of the background from that 
of the scatterer. The background contribution is described by the Bessel and Hankel 
functions, while the scatterer contribution is described by the function Xq- Standard 
multiple scattering homogenization theory is based on the asymptotic forms of all these 
functions to derive, by means of the monopolar and dipolar terms (q' = and g = 1), the 
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effective medium properties, as explained in [321 [33] • However, metamaterial behavior 
can appear, as it is demonstrated below, in the regime where only Bessel and Hankel 
functions of the background take their asymptotic forms. 

Let us then consider that the arguments of Bessel and Hankel functions are small 
[kbRa << 1), and use their asymptotic forms The monopolar component of the T 
matrix is 



2 """^o-'-a 2 

where the logarithmic term in the denominator is obviously equal to zero in the low 
frequency limit. However it cannot be neglected when dealing with metamaterials. This 
term has been omitted in most of preceding works about acoustic and electromagnetic 
metamaterials, but it contributes considerably to their effective parameters. 
Equivalently, the dipolar component of the T matrix is 

Since we expect this scatterer behaves as a homogeneous scatterer with acoustic 
parameters Pa and Ba, the matrix elements should have the form 



inRlkl 



(9 a) 



T, ^ 'ifP^L^kl (96) 

The comparison of these expressions with those given by equations ([7]) and ([8]) , permits 
to identify the frequency-dependent bulk modulus and density as 

Ba{uj)/B, = ^ In hRa - \hRaXo (10a) 

Pa{uj)/pb = Xl/hRa (106) 

Note that Ba{uj) and Pa{<^) are obtained from the fact that, after a scattering process 
(in the long wavelength limit) we expect to extract the parameters of a homogeneous 
fluid-like cylinder. These parameters are functions mass density at the surface of the 
scatterer, p{r = Ra), but they also depend on the field and its derivative at the surface, 
that is, of ipq{r = Ra) and dr'ipqir = Ra), respectively. These quantities are frequency- 
dependent and, consequently, they are the responsible of the frequency-dependence of 
parameters Baioj) and Paijjj)- 

In Appendix A we study three types of scatterers giving negative parameters at 
very narrow frequency regions. The first one is a homogenous fluid like cylinder such 
that Ca « Cb, condition that grants kb « ka- The second is also a homogeneous fluid- 
like cylinder but now with cylindrical anisotropy. These type of cylinders have already 
bee studied for cloaking devices, radial wave crystals [37] and hyperlenses [I6] and here 
we will see how they can be used to tune the resonance of the dynamical mass density 
of a metamaterial. Finally the third example shows that fluid-like shells can work as 
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Helmholtz resonators,obtaining from them negative bulk modulus, but also they can 
work as metamaterials with negative mass density. 

However, let us point out that we report here only the simpler examples of scatterers 
to realize acoustic metamaterials. Obviously more complex scatterers could be used to 
improve the frequency response, but a full analysis of this type is beyond the scope of 
the present work. 



4. Multiple Scattering of Acoustic Waves in the Quasi-Static Limit 

If we have a cluster of cyhndrical scatterers defined by some pa and Ba, we expect that, 
in the low frequency limit (that is, for wavelengths larger than the typical scatterer 
distance), they behave as a homogeneous medium with effective parameters pe// and 

The comparison among the scattering properties of the cluster and the effective 
scatterer is used to obtain the effective parameters for the case oflow filling fractions 
|43] . In [321 ESI [3l] such scattering formulation was generalized and all the multiple 
scattering interactions between the scatterers were included, dealing to more general 
expressions which also include the possibility of having anisotropy in the mass density 
of the effective medium. 

In the following subsections the results in [321 [331 El] are generalized to the case of 
metamaterials with frequency-dependent parameters. It is assumed that we can replace 
Pa and Ba by their corresponding frequency dependent values Pa{oj) and Ba{uj). It will 
be shown that this method is self-consistent and, therefore, it defines the correct way 
to explain acoustic metamaterials. 



4.1. Multipolar Interactions: The A Factor 

Let us consider that a cluster of scatterers are periodically arranged in the space. In 
the low frequency limit such a cluster behaves like an effective fluid-like medium with 
effective parameters given by 



5e//(u;) B^ Baioo) 

In the above expressions the dependence on the frequency has been added under the 
assumption that we can use for a scatterer the frequency-dependent parameters defined 
by equations fllOal) and (11061) . 

While equation flllap is valid for all filling fractions, equation (jllfel) is valid for 
diluted clusters only (i.e. low filling fractions). In [33j and [31] the expressions for 
the effective density were generalized for the case of high filling fractions, and all the 
multiple scattering terms were introduced in Equation (11161) by means of the so called 
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Figure 1. Frequency-dependent bulk modulus for a medium composed of air cylinders 
in a water background. 



A factor, leading to 

p,M(A + /)+MA-/) 
pa{uj){A - f)+pb{A + f) 

The A factor represents a multipolar correction to the effective density expression, 
and it includes all the multiple scattering terms between all the cylinders in a cluster or 
in an infinite lattice. Such a factor contains also information about the density of the 
cylinders forming the cluster, therefore if we want to define a frequency dependent 
A = A{(jj) factor and generalize the frequency dependent parameters to all filling 
fractions, the frequency-dependent mass density must also be included there. However 
this inclusion must be made carefully. 

However this factor is only important for very high filling fractions and also for 
very strong scatterers [32l [33], so we consider that such correction out of the scope of 
the present work; it only adds complexity to the calculation of effective parameters. 

4-2. Anisotropic Metamaterials 



Expressions in previous section have not taken into account the possibility of having 
anisotropy; that is, they did not consider the lattice symmetry in which cylinders are 
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Figure 2. Imaginary component of the effective speed of sound in a medimn composed 
of air cylinders in a water background. 



ordered in the cluster. In [31] expressions for the anisotropic mass density where found in 
terms of the cyhnder's parameters and the lattice geometry. It can be shown that, even 
neglecting the multiple scattering terms, we still have anisotropy for non-symmetric 
lattices, having the following expressions for the effective mass density tensor (the 
subindex "eff" has been omitted for clarity) 



where A is the anisotropy factor explained in Appendix and ri{u) is defined as 



It has been assumed that the lattice is oriented along the main axis and, therefore, the 
tensor has been previously diagonalized. 

The expression for the effective bulk modulus remains the same as in the previous 
section, so that we can obtain the effective sound speed tensor as [31] 



Note that the anisotropic mass density tensor can have both components with the 
same sign (negative or positive) or each component with a different sign. However, as the 
refractive index is the square root off those terms, if we want a negative refractive index 
we need a negative bulk modulus too, therefore we can have only or both components 




(13 a) 



(136) 




(14) 




(15) 
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Figure 3. Table summarizing the different propagation behavior of acoustic 
metamaterials. 



Figure 4. Effective parameters, normalized to those of the background, as a function 
of frequency of a rectangular lattice of fluid-like cylinders with b = 2a. The density, 
bulk modulus and radius of the cylinders are pa = O.bpb, Ba = 0.02Bf, and Ra = 0.49a, 
respectively. We see two narrow multipolar resonances between a/A = 0.3 and 0.35, 
however, in that limit the homogenization hypothesis is not good. The thin horizontal 
dotted line is a guide for the eye. 
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Figure 5. Real part of the effective speed of sound tensor of the medium described in 
Fig. 2] as a function of frequency. Note that we can have both components positive, or 
one negative and the other one imaginary (real part equal to zero in this plot). It is 
impossible have both components negative or one negative and the other one positive 
(see text). The horizontal dotted line is a guide for the eye. 



of the refractive index negative, or one negative and the other one imaginary (no wave 
propagation). In other words, we can have both an elhptical dispersion relation or a 
hyperbohc one. In Figure H] the anisotropic mass density and the scalar bulk modulus 
have been plotted, and in Figure [5] the corresponding components of the sound speed 
tensor are shown. 

5. Application to Electromagnetic Metamaterials 

The vectorial nature of electromagnetic (EM) waves makes the problem more complex, 
but in 2D the EM field can be decomposed into TE and TM modes, leading to the 
same wave equation as for scalar acoustic waves. Now P in equation ([1]) and ([5]) is the 
z component of the electric (magnetic) field for TE (TM) modes, and {p,B) = {fi,e~^) 
for TE modes and (p, B) = {e, fi~^) for TM modes. 

Although both problems are mathematically equivalent, physically they are quite 
different. The numerical values and ranges of the material parameters p, B and /x, e are 
not the same in both fields, therefore it is worth to study them apart each other. 

Thus, for TE modes 





In hRa + 



kaRa Jo{kaRa) £b 
2 Ji{kaRa)£a 



(16a) 



2 
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/Xf^(w) 1 Ji{kaRa)lJ^a 



fib kaRa J[{kaRa) fJ'b 

For TM modes, 



166) 



lnfc,i?. + h^MkaRa) lM, ^^^^^ 



£™H 1 MkaRa)ea 



(176) 



£b kaRa JlikaRa) 

These equations show that, even when the cyhnders are non-magnetic {fia = f^b = 
/io), we can have a magnetic response for a given frequency range. This is a well known 
phenomenon called mesomagnetism[H] 

Also these equations show that, as a function of frequency, the same cylinders 
behave in a different way for TE or TM modes, presenting different constitutive 
parameters that are equal when w — )■ 0, that is 

lim e^M = e^M = (18) 
lim/ir(u;)=/i™(c^) = /ia (19) 

The effective medium made of a cluster of these scatterers will present different 
constitutive parameters for each of the polarizations, so that, applying the results of 
Section HI 

ef//a;) = (1 - /)5, + (1 - /)e^^6 (20a) 
a^^ (u) - /fI!Mi±/)±i^^(i^^ (20b) 
/^™(^) = (1 - Z)/^^ + (1 - /)/^™& (20c) 

'^^f - .™M(i-/) + Mi + /)^^ ^^'^ 

Figure |6] depicts a plot of the effective constitutive parameters for a square lattice 
of dielectric cylinders with Sa = ll^b and Ra = 0.4a. Note that although f^a = fJ-b = A^o 
there is a strong magnetic resonance for both polarizations. However we see that the 
resonance of e-^*^ is beyond the homogenization limit (A < 4a), so that probably this 
resonance could not be observed. 

Figure [7] shows the effective speed of light (relative to that of the background) 
for this system. Note that only the TE polarization presents negative speed of light (or 
negative refraction index). This phenomenon is due to the fact that the resonance of 
e'^'^'^ is too far and too sharp to interact with that of /i^*^. If we still want to have 
negative refraction in the two polarizations we can use anisotropic cylinders, as we did 
for the acoustic case. 

If the cylinder's permittivity is given by a tensor of the form 

) (21) 
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Figure 6. Effective constitutive parameters, relative to those of the background, for a 
square lattice of dielectric cylinders such that Ra = 0.4a and = ll£b. The resonances 
of the permeability and permittivity for the TM case are so sharp that they do not 
allow to the system to present negative refraction for that polarization. However the 
TE polarization is allowed because of the wide range of negative permittivity. 



Figure 7. Effective speed of light normalized to that of the background for both TE 
and TM modes. As we can see, we cannot have negative refraction for the TM mode 
(see text). The regions where there is not solution for the speed of light corresponds 
to those in which e and /U have different sign, so that the speed of sound is imaginary. 
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Figure 8. Frequency shift to the left of e-^*^ due to anisotropy in the cyhnder. As 
we see, there is a frequency region in which ah the constitutive parameters for the two 
modes (TE and TM) are negative, so that we have a region of total negative refraction. 



the expressions for the frequency-dependent constitutive parameters are now 



(226) 



-*M,„,A + ^J!!^^4™^^i (22c) 



(22 



2 " " ' 2 Ji(fc™i?.)/i. 

en^) _ 1 Mk™Ra) eae 
e, k™Ra J^,{k™Ra) e, 

where /cj^ = Uy/e^^, fcj*^ = u^e^a and 7^ = Eae/eav 

These expressions allow us to shift the resonance of e^^^ to lower values just by 
increasing the e^r, keeping the rest of the system unaltered, in the same way as we did 
for the acoustic case. Note that the equivalence is pr — )■ ^^^^ and pe — )■ e^*^, that is we 
need to increase the value of Ear in the electromagnetic case. Thus, as we increase it, 
the anisotropy factor 7 goes to 0, and the resonance moves to the left. 

Figure E] consideres the same system as in figure |6] but with anisotropic cylinders, 
where Sr = 200eg and eg = = llSb- Note how all the resonances keep their positions 
but e™, which now moves to the left reaching , leading therefore to an effective 
medium with negative refraction. The effective speed of light is depicted in Fig[9l where 
it is obvious now that both polarizations present negative refraction properties within 
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Figure 9. Effective speed of light for both TE and TM modes for the system of Fig. 
[51 It is clear that there is a region of total negative refraction, thanks to the frequency 
shift of e-^*^ due to anisotropy (see text). 

a similar frequency range. 
6. Summary 

In summary, a comprehensive multiple scattering formulation of acoustic metamaterials 
has been introduced. This formulation is based on a homogenization theory in the 
quasi-static limit, in which we allow the wave number in the background be arbitrarily 
small while the wave number in the scatterers remains finite. 

In general, it is shown that ordered or disordered arrays of sound scatterers can 
behave, in the low frequency limit, like effective fluid-like materials with either positive 
or negative acoustic parameters, where these acoustic parameters are the scalar bulk 
modulus and the tensorial mass density. The behavior of the effective medium depends, 
among other properties, on the surface fields in the scatterer. Therefore, it is possible 
to improve or manage the frequency response of the effective medium with complex 
scatterers, like fluid-like shells or anisotropic fluid-like materials. 

Examples of these scatterers have been analyzed, showing that they present 
negative effective parameters whenever the theory predicts them, verifying therefore 
the formulation presented. Also, it has been shown how these complex scatterers can 
be used to tune the effective parameters of the medium. 

The theory developed for acoustic waves has been extended to electromagnetic 
waves in 2D, showing that equivalent type of scatterers can also be used in 
electromagnetism to tune the effective medium response. 
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Figure Al. Frequency-dependent acoustic parameters of a homogeneous cylinder 
with Ba = 0.005-Bb and pa = 0.5p(,. The cyhnder's radius is Ra = 0.3a, where a is the 
lattice constant. 

In conclusion, the theory presented not only explains the metamaterial behavior 
found so far in the literature, but also gives the basis for improve this behavior with 
more complex scatterers. 
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Appendix A. Examples of Scatterers with Local Negative Parameters 

Appendix A.l. Homogeneous and Isotropic Scatterers 

For a homogeneous scatterer with parameters Pa and Ba the field inside the scatterer 
is given in terms of Bessel functions. Therefore, after some algebra, the frequency 
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where /ca = uj^'p^jB'a. For fc^ — we recover the cyhnder's parameters. Note that the 
metamaterial concept appears when k\j «\ but is not small enough. This condition 
occurs when Ca « c^; i.e. when the wavelength in the background is several times 
larger than that inside the scatterer. 

Figure lAll plots the frequency-dependent parameters described by Eqs. f llOap and 
( ]106[) for a soft scatterer with Ba = 0.005-B;, and pa = O.Spf,. These values give a sound 
speed of the scatterer Cq = 0.1c{,, which locates the resonances in the low frequency 
limit, as can be seen in the figure. 

From (1106P it is deduced that the region of negativity for Pa{oj) is determined 



by the first zeros of Ji{kaRa) and J[{kaRa), which are a = 3.8317 and a' = 1.8412, 

respectively |36j, and corresponds to frequencies co- such that 

1.8412c, ^ ^ 3.8317c„ . 2c, ,^ 
B < uj^ < , Aw_ ^ — (1.2) 

-^a 

Therefore, if we want to locate the frequency region in the low frequency limit, we have 
to decrease c, (for fixed Ra). But, as a consequence, the bandwidth will decrease; that 
is, the resonance becomes sharper. 

These type of scatterers are possible only in a dense background, like water, where 
we can get such a low bulk modulus and density. If we need to get metamaterials in an 
air or gas background another approach should be used instead. 



Appendix A. 2. Homogeneous and Anisotropic Scatterers 

Sometimes it is not possible to obtain sound speed smaller than a certain value but we 
still want to decrease the frequency at which negative mass density appears. In this 
case, fluid-like cylinders with circular anisotropy |38] can be used to shift the resonance 
to lower frequencies. 

Anisotropic cylinders are characterized by a scalar bulk modulus Ba and a tensorial 
mass density whose componentes are constant when referred to a cylindrical coordinate 
system, = (p^, Pe)- In these cyhnders the pressure field is described in terms of Bessel 
functions of real order 7g, where 7 = \fihTpe is the anisotropy factor, 

^g{r,u}) = Jyg{kar) (1.3) 

Despite being anisotropic cylinders, the circular symmetry of these scatterers make 
them suitable for applying the theory developed in this work. Thus, when q = the 
field distribution is the same than that of an isotropic cylinder (because J^q for g = 
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Figure A2. Frequency-dependent dynamical mass density of an anisotropic fluid like 
cylinder for three different anisotropy ratios. The radial component of the sound speed 
tensor is Cr = 0.2c{, in the three examples. We see that, as we increase the anisotropy 
ratio, the resonance is shifted to lower frequencies. 



is Jq); therefore the frequency dependent bulk modulus will be also given by (11. lap . 
However ( 11.161) now becomes 



and, as we let the anisotropy factor 7 be smaller, the dipolar resonance g = 1 becomes 
closer to the monopolar one. 



decreasing therefore the resonant frequency. 

In Fig. IA2l the frequency- dependent mass density is plotted for three different Pr/ Pe 
ratios. Note how as we increase the value of pe (so that we decrease the anisotropy factor 
7) the resonance of the density moves down in frequency. 

These type of strongly anisotropic fluid-like cylinders have already been used by 
Li et al. [I6] for building acoustic hyperlens and it has been recently characterized for 
different anisotropy ratios in [39], showing also the same frequency shift. 

Appendix A. 3. Fluid-Like Shells as Helmholtz Resonators 

Let us assume now that we have a fluid like cylinder of radius Ra and parameters pa 
and Ba- If this cylinder is enclosed by another of radius Rh > Ra and parameters ps and 
Bs we have a fluid-like shell. Obviously this is an idealization, because such a structure 





7^0 
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Figure A3. Frequency-dependent bulk modulus of a water shell in air with ps = 
lOOOpb and Cg = 4.4cf,. The outer radius of the shell is Rt ~ 0.3a, and the plot shows 
three different radii Ra- 



cannot be realized with common fluids. However, on the one hand, if the shell is an 
elastic material, this can be sometimes a good approximation and, on the other hand, 
if the fluids are "metafluids" [iQl ST], |12] , the structure can be easily fabricated. 

These structures can work, as Helmholz resonators [2H EE]. Here we give a more 
rigorous derivation of the resonance frequency. 

After applying the correct boundary conditions, the impedance factor Xq of a fluid- 
like shell is given by 

p,cs UkA)+T-Y,{KR,) 

Xq = - 



where 



and 
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The parameters of the shell are ps and c^, while those of the enclosed cylinder are pa 
and Ca- 

For g = and fc^ — ^ the inner T matrix Tq is equal to 
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where 



Ba 2 

Thus the impedance factor of the shell can be approximated to 
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which can be zero only for Bg/Ba » 1 and consequently 
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which defines a cut off frequency 

Cn 



OJr. = 
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If the shell is soft, that is, if Bg/Ba << 1, a negative bulk modulus appears as in 
the homogeneous cylinder case, and the shell nature of the scatterer is not relevant. In 
that Rb — > Ra the impedance factor reduces to 
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(1.14) 



and the bulk modulus cannot be negative, however we will see now that, in that case, 
the density can be negative. 

When q — 1 the inner T matrix is approximated by 

T^Rlks Pa - Ps 



rpa , 

^ ' 4 Pa + ps 

the density becomes negative once the denominator of xi cancels, that is, when 



J[{ksRb)+T^Yi{ksRb 



k;Rb Ra Pa - Ps 



4 Rlpa + ps 

This expression gives a cut off frequency for the negative density of 



= 0. 
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where now we need that ps « Pa- 



Appendix B. Technical Details 
Appendix B.l. The anisotropy factor A 

A two dimensional periodic array of scatterers is defined by the lattice vectors ai and 
a2, so that the position Rn of any scatterer in the lattice can be determined by two 
integers rii and 77-2 such that 

Rn = niai + n2a2 (2.1) 
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This lattice has also associated the reciprocal lattice vectors 61 and 62 such that 

b, ■ aj = 27r6ij , i = l,2 (2.2) 

If we define the reciprocal lattcie point Gh = {Gh, Oh) as 

Gh = hbi + h^b^ (2.3) 

the anisotropy factor A can be found in [31] and is given by 

A = AS :Mg^!!^e-2^^'^ (2.4) 

where Jz{') is the third order Bessel function and Rmin is the smaller of bi, 62? ^1 + ^2- 
Factor A can be made always real by properly choosing a coordinate system in which 
the tensors be diagonal. 
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